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(54) Devices for in-vivo determination of the compliance function and the systemic blood flow of 
a living being 



(57) A device for individual in-vivo determination of 
the compliance function C(p) = dV/dp of the vascular 
system downstream of a ventricle and/or systemic blood 
flow of a living being from the Wood pressure p(t) and a 
reference cardiac output COref. The device comprises 

a) a pressure sensor which continuously deter- 
mines the pressure p(t) in the aorta or in the vicinity 
of the aorta; 

b) a computer which 

bl) calculates the mean blood pressure MAP 
from the blood pressure p(t), 
b2) calculates the systemic resistance R of the 
body as 



b4) calculates the compliance function C(p) at 
least from p(t), p(t) and R using a nonlinear 
model. 



R = 



MAP - CVP 
COref 



CVP being an arbitrary central venous pres- 
sure which is ascertained or estimated, and 
COref being a reference value for the cardiac 
output, 

b3) takes at least the first differential of the 
blood pressure with respect to time 
p(t) = dp/dt,and 
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Description 

[0001] The invention relates to a device for individual in-vivo determination of the compliance function C(p) = dV/dp 
of the vascular system downstream of a ventricle of a living being from the blood pressure p(t) and a reference cardiac 
output COref . 

[0002] The invention also relates to a device for continuously determining the systemic blood flow q(t), in which the 
blood pressure p(t) in the aorta or in the vicinity of the aorta is determined continuously. 

[0003] Methods and devices of the type mentioned above are known. In the past, a plurality of methods have been 
developed with the purpose of calculating the cardiac output (CO) from the arterial blood pressure. On the one hand, 
methods are proposed in which the CO is determined from a few characteristic values such as, for example, mean arte- 
rial pressure (MAP), systolic and diastolic pressure (APsys, APdia), ejection time (ET) and patient data (age, sex, 
weight, height) [4,6,7], and on the other hand algorithms are used in which the complete contour of the pulsating blood 
pressure curve is utilized to calculate the cardiac output [1,5,20]. In the latter methods, which are also referred to as 
pulse contour analysis, two different approaches have so far been adopted. Firstly, the CO has been calculated directly 
from the arterial blood pressure with the aid of some correction factors [19,20] while in other work [5,25] a blood flow is 
calculated from the pressure, together with particular assumptions, and is assumed to correspond to the actual blood 
flow in the aorta and therefore to be usable for determining the cardiac output. 

[0004] The classical Windkessel model, which was first proposed by Hales [26] and has been used by Frank [27] to 
determine the stroke volume (SV) and, together with the heart rate, the cardiac output, uses only the peripheral resist- 
ance R and the compliance C for modeling the cardiovascular system in question. In this model, the arterial blood flow 
is described by q(t), which can be calculated for given C and R with the aid of the blood pressure p(t) which is to be 
measured. However, closer examination shows that this simple model reproduces the physiological conditions only 
incompletely, with the result that many modifications to the original model have been proposed; for an overview refer- 
ence may be made to [24,28]. 

[0005] The accuracy of these methods depends essentially on how well the assumptions, i.e. the model used, reflect 
the conditions in the cardiovascular system in question, and in [5] a nonlinear Windkessel model is thus used whose 
parameters are dictated by the age and sex of the patient. More recent investigations [21] show, however, that the 
model used in [5] does not reproduce the correct underlying physiological conditions; in particular the compliance 
(extensibility) of the vessels cannot always be described by the compliance/pressure relationship given in [5]. There 
may be several causes for this discrepancy. First, only a dependence of the in-vitro determined aortic cross section on 
the blood pressure is taken into account in [5] and a length variation, as described in [22,23] is neglected; also, only the 
density of the blood and not the strongly hematocrit-dependent viscosity is taken into account, and the compliance of 
the peripheral system is likewise ignored. Secondly, apart from age and sex, the compliance function C(p) of a particu- 
lar individual cannot be used in the method described in [5]. However, it is precisely in the examination of pathophysio- 
logical cases, e.g. arteriosclerosis, that it cannot be assumed that C(p) varies according to age and sex, so that the 
basic model describes the physiological conditions only incompletely [25]. Lastly, it has been shown in [24] that it is to 
be expected that a modified Windkessel model can reproduce underlying physiological conditions more precisely. 
[0006] However, a common factor in all the models described above is that the model parameters, after they have 
been determined once, no longer depend on the condition of the cardiovascular system in question. Nevertheless, 
almost all parameters can change wfth time, and for example the systemic resistance R can change as a result of med- 
ication. Other parameters, including the expandability and length of the aorta, change so greatly with pressure that they 
actually have to be regarded as variable even within one heartbeat. 

[0007] The fact that aortic impedance and compliance cannot be assumed to be constant has been shown both in 
animal experiments [22] and for humans [29]. Primary causes of this are the expandability, length variation and volume 
variation of the aorta and vessels in proximity to the aorta. The typical variation in the aortic volume V as a function of 
pressure has been described inter alia in [30]. Since the compliance of the system is given by 



and because of the limited volume the compliance must tend toward zero for very high pressures and cannot be con- 
stant. Since the change in volume is caused by length and cross-sectional changes in the vessels, there is also a 
change in the aortic impedance which, according to the Navier-Stokes equation, is determined on the one hand by the 
cross section and density of the liquid and, on the other hand, by length, viscosity and density of the blood. 
[0008] Pressure-dependent aortic impedance and compliance have been discussed inter alia in [5,21] and used 
therein to investigate nonlinear Windkessel models. In [5] it is in particular assumed that C(p) can be established by age 
and sex of the patient. The impedance Z(p) is also ignored in this approach. What is more, it follows from the results 
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obtained in [21] that the model used in [5] may to some extent conflict with the true physiological situation. One cause 
of this is that the compliance and aortic impedance are preset. An approach of this type is unsuitable for taking into 
account the features characteristic of the patient in question. In addition, the method proposed in [5] cannot be applied 
without modifications to other species. Further, only the typical aortic diameter investigated beforehand in-vitro and the 
s density of the Wood are taken into account in [5J. The effect of aortic length variations, and the dynamic behavior of the 
vessels in proximity to the aorta and the peripheral vessels and the viscosity of the blood are ignored in the modeling 
of the conditions existing in vivo. 

[0009] There is consequently not yet any method which, for individual in-vivo determination of the compliance/pres- 
sure relationship, employs the measured variables used here. 

10 [0010] These disadvantages are to be eliminated by the device of the present invention by determining all model 
parameters of interest from measurements on the physiological system in question, i.e. human or animal. To this end, 
in particular, the blood pressure p(t) in the aorta or in proximity to the aorta is to be measured continuously and a refer- 
ence cardiac output (COref) is to be measured at least one time. With the aid of these values, all the parameters can 
be established and then used for hemodynamic investigation. 

is [001 1 ] The object of the invention is to provide a device for individual in-vivo determination of the compliance function 
of a living being, which reproduces physiological conditions as faithfully as possible. 

[0012] A further object of the invention is to provide a device for continuously determining the systemic blood flow of 
a patient which has a low degree of invasiveness and describes the actual blood flow at any given time as accurately 
as possible. 

20 [0013] The way in which these objects are achieved is described in the independent patent claims. Advantageous 
refinements can be found in the dependent patent claims. 

[0014] According to the invention, a so-called Windkessel model is used whose parameters can be identified with the 
aid of a reference measurement in vivo. Subsequently, the systemic flow and other hemodynamic parameters are 
thereby determined. A Windkessel model adapted and modified in this way describes the cardiovascular system of the 

25 individual in question more accurately and can therefore be used for likewise more accurate calculation of the systemic 
flow and hemodynamic parameters derived therefrom. The method can also be applied directly to other species, with- 
out the need to determine previously hypothetical characteristics for this. The extra outlay for the newly developed 
method consists in the fact that, for calculating an individual compliance function, besides the continuous blood pres- 
sure measurement, the cardiac output has at least one time to be determined using a different method, e.g. arterial 

30 thermodilution. 

[0015] The invention will be explained in more detail below with reference to an illustrative embodiment schematically 
represented in the figures, in which: 



Figure 1 shows a preferred electrical model circuit for simulating the cardiovascular system in question. 
35 Figure 2 shows a flow chart for calculating the aortic impedance Zip), the compliance function C(p) and the blood 
flowq(t). 

Figure 3a shows a graph to illustrate the time dependence of the blood pressure p(t), to denoting the time at which 

the aortic valve opens, ts the time at which the aortic valve shuts and tD the end time of the diastole. 
Figure 3b shows a graph to explain the time dependence of the blood pressure p(t) and the resulting blood flow q(t). 
40 Figure 4 shows a typical compliance function C(p) of a human aorta. 

Figure 5 shows a block circuit diagram of a device according to the invention. 



[0016] Figure 1 shows a nonlinear modified Windkessel model which is preferably used according to the invention, in 
which the aortic impedance functions Z(p) and Zo(p), the compliance function C(p) and the systemic resistance R are 
45 taken into account. 

[0017] The resistor R in Figure 1 represents the time-variable peripheral flow resistance of the body. Zo(p) and Z(p) 
are nonlinear impedances which are dependent on the pressure p(t) and, together with the nonlinear pressure-depend- 
ent capacitance G(p), are intended to simulate the behavior of the aorta and the blood vessels in proximity to the aorta. 
[0018] The result obtained for the model outlined in Figure 1 with the Fourier transform 

50 

00 

p(a>) = J pWextf-ifflQdt 

-ao 

55 

and the function q (©), which is to be calculated similarly, for the aortic impedance is 
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. . q(co) 11 

lira —— = - + - , (1) 

<0-k*> p (fl)) K i> 



so that for R > > Z it follows that 

10 



p«fl) 

Z - lim z — • (2) 

«-♦- > q(CO) 



75 

[001 9] The following equation is further satisfied for the compliance C 



C (n\ - Q(t)-P(t)/R n x 
p(t)-2(p).(q(t)-p(t)/R)' w 



20 



in which p(t) = dp(t)/dt and q(t) = dq(t)/dt are the respective differentials of the pressure and flow with respect to time. 
Equations (1) to (3) show that C and Z can be calculated for the model in Figure 1 if the systemic flow q(t), the blood 
pressure p(t) and the systemic resistance R are known. However, the model describes the cardiovascular system in 
25 question well only to the extent that this is possible in the scope of the proposed approach. 
[0020] Figure 2 shows an overview of the refined method. 

i) The mean arterial Wood pressure MAP and the heart rate HR are firstly determined from the pressure p(t). 

ii) Together with the reference cardiac output COref, which has preferably been determined by arterial thermodilu- 
30 tion, and for which the following equation is satisfied 

COref = HR-Jq(t)dt (5) 

the systemic resistance is calculated according to R=(MAP-CVP)/COref . In this expression, CVP is the central 
35 venous pressure which, if ft is unknown, can be approximated by a suitable constant pressure, e.g. 0 mmHg. 

iii) The next step is to establish a flow q(t), which should be chosen suitably, which is used as the start function in 
the subsequent iteration and satisfies the underlying physiological conditions. The blood flow q(t) describes the flow 
which passes directly from the left heart into the aorta. It is therefore required of q(t) that the subconditions 

q(tO)=q(ts) = 0 and Jq(t)dt =* ^ = SV (6) 

are satisfied for all heartbeats recorded during the measurement of the reference cardiac output COref. In this 
45 expression, the time to is the start of the systole and ts the end of the systole. The end of the diastole is denoted 
tD below. Suitable flows q(t) would be q(t) = 0 for ts < t <> t D and for to < t ^ ts, for example 



or preferably 
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q(t) - KtP(t)-(p(to)*(t-to) P(t, t |'.^ (tQ} ) I . (8) 

5 

in the latter case, k is to be determined from the condition (q(t)dt = COref/HR . Besides the flows indicated in equa- 
tions (7) and (8), however, other starting conditions are conceivable, for example a constant or parabolic flow. 

iv) Set Z(p) = 0 and introduce the auxiliary variables qold(t) and Eold, which are initialized at qold(t)=q(t) and 
10 Eold = oo . 

v) Calculate, from the blood pressure p(t) measured in the aorta or in proximity to the aorta, the flow q(t) and its 
time derivatives, the compliance function according to 



15 



25 



35 



40 



Cfpi. <ffl-PW* 

p(t).«p)Mq(t).p(t)/H) . ( 9 ) 



20 vi) The inverse of the compliance function is then approximated by a polynomial of suitable order in p, i.e. in 

i ■ Z a ^ k < 10 > 



CM k 



the expansion coefficients are determined in such a way that the identity 

q(t) - 3® + • 1 k -- m -VP) '(q(t) • pdVRH 

* C£fc«fcP*] (11). 



is optimally satisfied. As a suitable criterion for this, the mean square error 

E a ((q«.pttffl'|&C).Z(tf .(*!)■ pWR)l /BkOW^l)" ) 



may be minimized, it being possible to use p(t) and q(t) at all times or alternatively only from preferred time intervals. 
In this case and below, the notation < • ) indicates calculation of the mean. 

vii) If E < fold, then set qold(t) = q(t) and Eold = E and continue with step viii), otherwise go to point x). 



45 viii) Calculate Z(p). On the one hand, the procedure adopted for this may be to determine Z(p) according to equation 

(1) or (2). In this case, it is assumed that these equations, which are initially valid only for the model in Fig. 1b in which 
the parameters are not pressure-dependent, also apply for sufficiently short time intervals At to the nonlinear approach 
proposed according to Figure 1 . For the latter, the impedance function can then be ascertained according to equation 

(2) with 



♦At 

Z(p= Jp(t)dt/2At) 

55 or, equally, by 
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10 



J t+At t+At t+At t+At 

( J p 2 (t)dM J p(t)dt) 2 )/( J q 2 (t)dt-( J q(t)dt) 2 ) (12) 
t-At t-At t-AI t-At 

even directly from the time-dependent blood pressure p(t) and the Wood flow q(t) without a prior Fourier transform. 
On the other hand, Z(p) can be calculated yet more simply by 



as can be seen from the following discussion. The aortic diameter d and length I increase with rising pressure, so 
that to first approximation d I may be assumed. According to the Hagen-Poiseuille law, this results in Z(p) a VV, in 
is which t| denotes the viscosity of the Wood and V the aortic volume. With C(p) = dV/dp , this gives C(p) a d(1/Z)/dp, 
from which equation (13) directly follows. The constant of proportionality A which it contains can, for example, be 

determined by determining the function Z(p) for at least one pressure p according to equation (1 2). 

ix) If the error E is small enough, then the identification of the model parameters is terminated here, Otherwise, it 
continues with x). 

20 x) The assumed blood flow should finally be varied in such a way that the stroke volume further corresponds to the 
stroke volume SV = COref/HR which follows from the reference cardiac output. Since, at this time, qold(t) always 
describes the optimum flow so far, q(t)=qold(t) + 8q(t) with J5q(t)dt=0 is set. 
xi) Continue with step v). 



25 [0021 ] The algorithm indicated in i)-xi) describes the preferred method, in which the reference cardiac output COref 
and the continuously measured arterial pulse curve p(t) are used to determine all the other values. This guarantees that 
the compliance function and the aortic impedance function are determined in such a way as is required by the interac- 
tions actually taking place in the cardiovascular system in question. In particular, C(p) thereby takes into account not 
only the variation in the aortic cross section, but also the actual variation in volume of the aorta and peripheral vessels; 

30 likewise, a variation in the length of the aorta and the density and viscosity of the blood are described by Z(p). With the 
aid of the description which is introduced in step vi) and can reproduce any physiologically feasible compliance/pres- 
sure relationship, it also becomes possiWe to extrapolate C(p) and Z(p), so that the functions can even be applied 
beyond the pressures observed during the calibration, i.e. during the measurement of the reference cardiac output. 
[0022] The algorithm indicated describes the preferred method. Other methods can be readily obtained from it and 

35 are covered by the appended claims. In particular, in vi), instead of the inverse of the compliance function, it is also pos- 
sible for C(p) to be described by a finite Taylor series, i.e. by a polynomial. It is also conceivable for steps viti) and ix) to 
be interchanged with each other, or for the criterion to be optimized to be modified in step vi). For example, instead of 
minimizing the mean square error, it is also possible for the expectation to be maximized. 

[0023] In order to accelerate the method, it is in particular possible to configure the choice of the initial Wood flow q(t) 
40 in step iii) in such a way that C(p) is initially determined only from the diastoli, i.e. for ts < t <, Xq. According to equation 
(9), C(p) = - p(t) / R • p(t) is then satisfied for Z=0, so that for all p(t) < p(ts) the compliance function can therefore be 
used to expand the flow during the systole in terms of an orthogonal function system, e.g. 

45 q«-Ifl^^-U%y) , (14) 



The coefficients qk are in this case determined by minimizing the expression 

so <((q(t)^)(t)/R-p(t)/C(p)) 2 > (15) 

the above error being to be calculated from all p(t) < p(ts), including those which occur at the start of the systole. With 
this choice of the flow, if the cardiovascular system is in good enough condition, just one or a few iterations will be 
enough to calculate the model parameters using the described algorithm. 
55 [0024] The method may be further accelerated in the case in which, in step vi), only pressure values from the diastole 
ranges are used and 1/C(p) is approximated by a second-degree polynomial. During the first iteration, the following 
result is obtained for this application 
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cr . MAP** COff 1 

<P(t)> 3'MAP 'P-3 'MAP a .p 2 (16 ) 



[0025] After the parameters of the model in Figure 1 have been determined from the blood pressure p(t) in the aorta 
or in proximity to the aorta and from the reference cardiac output COref, the cardiac output can subsequently be calcu- 

10 lated continuously just from the blood pressure beat by beat. For this purpose, the start of the heartbeat, which is again 
denoted to below, needs to be calculated and equation (11) needs to be calculated for the initial condition 
q(t 0 ) = Oft o) =0 • A calculation of this type can be carried out numerically e.g. using Runge-Kutta, Euler or another gen- 
erally known algorithm. If, further, the required blood pressure measurement is available only at discrete time intervals, 
then linearization can be carried out between them if required. For the stroke volume thus determined only using the 

is blood pressure, the following is satisfied with y = 1 

SV = yJq(t)dt (17) 

so that the following is obtained for the continuously determined cardiac output CO 

CO = HR • SV = HR • y • fq(t)dt (1 8) 



20 



[0026] In the expression above, HR denotes the heart rate, which is likewise to be calculated from the pulse curve 
p(t) and indicates the number of beats per minute. The integration occurring in equations (16) and (1 7) can in this case 

25 be carried out over the entire heartbeat, or alternatively only over the length of the systole, since q(t)=0 is satisfied dur- 
ing the diastole. If the stroke volume SV, and therefore CO as well are calculated during the entire heartbeat, then it is 
not necessary to determine the end of the systole. To do this (see e.g. [31]), either accurate analysis of the pressure 
curve would be necessary, in order to determine the position of the so-called dichrotic notch and therefore the end of 
the systole from p(t), or further measuring instruments such as an ECG would be required. Integration over the entire 

30 period is consequently more robust and less involved than those methods which evaluate only a particular period of the 
heartbeat. If, in addition, the continuously determined cardiac output CO is also calculated from those blood pressure 
measurements which, together with the reference cardiac volume COref, have been used for model identification using 
the method described above, then the accuracy of the method for continuous CO calculation can be increased further 
in that CO=COref must be satisfied and the calibration factor y is therefore determined according to 

35 

COref . . 

Y= HR-jq(t)dt 

40 [0027] In order to use the described method, it is necessary to have a device whose basic structure is represented in 
Figure 5. In this diagram, the components represented by a broken line are optional and at least some of them may be 
omitted in a minimum configuration of the device. A device of this type consists at least of an evaluation unit, usually a 
central processing unit, in which the method for determining the individual compliance function C(p) is implemented on 
its own or together with other methods; in particular, the method for continuously calculating the cardiac output may be 

45 employed in the same device. It is also necessary to have a sensor for measuring the blood pressure p(t) in the aorta 
or in proximity to the aorta and an arrangement for signal processing and signal conversion, a program memory and a 
data memory, as well as a device for providing the reference cardiac output COref. If COref is determined through arte- 
rial thermodilution, then this unit consists at least of a blood temperature sensor and a sensor for measuring the tem- 
perature of the injection dose used by this method, see [8]. Since, however, COref can also be obtained in other ways, 

so this module may also have a different form or entry may take place through a keypad which may also be used in the 
device for the user to enter instructions. There will also be at least one of the options to have the results of the evaluation 
displayed, printed out or stored on a bulk storage device (not shown). 



55 
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Claims 

1 . A device for individual in-vivo determination of the compliance function C(p) = dV/dp of the vascular system down- 
stream of a ventricle of a living being from the blood pressure p(t) and a reference cardiac output COref, wherein 

30 the device has the following : 

a) a pressure sensor which continuously determines the pressure p(t) in the aorta or in the vicinity of the aorta; 

b) a computer which 

35 b1) calculates the mean blood pressure MAP from the blood pressure p(t), 

b2) calculates the systemic resistance R of the body as 

R MAP - CVP 
COref ' 

40 

CVP being an arbitrary central venous pressure which is ascertained or estimated, and COref being a ref- 
erence value for the cardiac output, 

b3) takes at least the first differential of the blood pressure with respect to time p(t) = dp/tit , and 
45 b4) calculates the compliance function C(p) at least from p(t), p(t) and R using a nonlinear model. 

2. The device as claimed in claim 1 , wherein the computer only uses values of p(t) which satisfy the following condi- 
tion for calculating the compliance function C(p): 

6 P(t)*P(tS), 

ts being the time when the aortic valve shuts. 

3. The device as claimed in claim 1 , wherein the computer only uses blood pressure values from the diastole for cal- 
55 culating the compliance function C(p). 

4. TTie device as claimed in claim 1 , wherein the computer only uses blood pressure values from the systole for cal- 
culating the compliance function C(p). 
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5. The device as claimed in one of the preceding claims, wherein the computer determines a blood flow q(t) on the 
basis of the pressure p(t) and the first time derivative dp/dt, and calculates the compliance function according to 

5 m p(t)-Z(p).(q(t)-p(t)/R) 

for arbitrary impedance functions Z(p) and arbitrary times t in such a way that 
w Off) = $ + C(p)[p(t) - 2(p) • (q(t) - p(t)/R)] 

is optimally satisfied. 

6. The device as claimed in one of the preceding claims, wherein the computer approximates the inverse of the com- 
15 pliance function, i.e. 1/C(p), by a finite-order polynomial and uses this polynomial to extrapolate C(p) beyond the 

pressure range recorded when determining the reference cardiac output. 

7. The device as claimed in claim 1 , wherein the computer deteTmines the minimum of the function 



20 



30 



45 



50 



z V f q(t) . m . p(t)-z(p)-(g(t)-p(t)/R)i 2dt 

t U R S k a k p k J 



25 and then calculates the individual compliance function C(p) as 

1 



C(p)= 



I k a kP k 



8. The device as claimed in claims 3 and 6, wherein the computer describes the inverse of the compliance function 
C(p) by a second-order polynomial and approximates C(p) by the following function: 

c(p)= MAp2c0ref - 1 



(p(t)> 3* MAP -p-3 • MAP 2 -p 2 

9. The device as claimed in one of the preceding claims, wherein the computer approximates the compliance function 
C(p) by a finite-order polynomial and uses this polynomial to extrapolate C(p) beyond the pressure range recorded 

40 when determining the reference cardiac output. 

10. The device as claimed in claim 1 , wherein the computer determines the minimum of the function 

tot, 

£ J (q(t) - ^ - [I k p k p k ] • [p(t)-Z(p) • (q(t)-p(t)/R]) 2 dt 
t-t e 



and then calculates the individual compliance function C(p) as 

C(p) = E k p k p k 



11. The device as claimed in claims 3 and 5, wherein the computer uses the compliance function C(p) calculated for 
p(t) <, p(ts) to expand the blood flow q(t) in terms of a complete function system and, in particular, describes q(t) in 
55 the form of a Fourier series by the following equation 
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q(t) = £q k sin(k-«. t -^. 

5 the coefficients qk being determined by minimizing the mean square error and the values to and ts denoting the 
times when the aortic valve is opened and shut. 

12. The device as claimed in one of the preceding claims, wherein the computer varies the assumed blood flow q(t) in 
such a way as to minimize the mean square error. 

10 

13. The device as claimed in one of the preceding claims, wherein the computer determines the aortic imped- 
ance/pressure relationship by 



15 




A being a constant of proportionality. 

20 14. The device as claimed in one of the preceding claims, wherein the computer determines a nonlinear aortic imped- 
ance function using the Fourier transforms of the blood pressure p {<£>) and postulated blood flow q (©) according 
to 1/Z = q(©->ao)/p(co->oo)-1/R or Z(p) = p(oo->ao)/qco->°o) 

15. The device as claimed in claim 1 , wherein the computer determines a blood flow q(t) on the basis of the pressure 
25 p(t) and the first time derivative dp/dt and calculates the impedance function 

+At 

Z(p= J p(t)dt/2At) 

-At 

30 

according to 

__ _ _ _ 

35 Z(p)= ( I p 2 (t)dt-( J p(t)dt) 2 )/( J q 2 (t)dt-( J q(t)dt) 2 ) 

^ t-At t-At t-At t-At 

16. The device as claimed in one of the preceding claims, wherein the computer approximates the aortic impedance 
40 function 2(p) by a finite-order polynomial and uses this polynomial to extrapolate Z(p) beyond the pressure range 

recorded during calibration. 

17. A device for continuously determining the systemic blood flow q(t) of a living being, wherein the device has the fol- 
lowing: 

45 

a) a pressure sensor which continuously determines the pressure p(t) in the aorta or in the vicinity of the aorta; 

b) a computer which 

b1) calculates the mean blood pressure MAP from the blood pressure p(t), 
so b2) calculates the systematic resistance R of the body as 

D MAP • CVP 
H= COref ' 

55 

CVP being an arbitrary central venous pressure which is ascertained or estimated, and COref being a ref- 
erence value for the cardiac output, 

b3) takes at least the first differential of the blood pressure with respect to time p(t) = dp/dt , and 
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b4) calculates the compliance function C(p) and the blood flow q(t) at least from p(t), p(t) and R using a 
nonlinear model. 

18. The device as claimed in claim 1 7, wherein the computer determines the systemic Wood flow q(t) by: 

5 

q(t) = ejS + C(p)[p(t)-Z(p) . (q(t)-p(t)/R)] . 



io 1 9. The device as claimed in claim 1 7, wherein the computer uses the compliance function C(p) ascertained according 
to one of claims 1 to 16. 

20. The device as claimed in claim 1 7, wherein the computer calculates the stroke volume SV by integrating the flow 
over a suitable period of time in accordance with 

15 

SV =/q(t)dt 

it being in particular possible for the suitable period of time to correspond to the heartbeat or the ejection time dur- 
ing the heartbeat. 

20 

21. The device as claimed in claims 1 and 17, wherein the computer calculates the stroke volume SV, by comparing 
the continuous blood flow q(t) with a reference cardiac output COref , by 

SV = y f Q («)dtwith y = R ^g s 



22. The device as claimed in one of claims 16 and 17, wherein the computer calculates the stroke volume variation 
according to 

30 

sw=V<sv 2 >-<sv> 2 

and uses it on its own, or with other parameters, e.g. mean blood pressure MAP, systolic pressure APsys, diastolic 
pressure APdia and heart rate HR, to correct the stroke volume. 

The device as claimed in one of the preceding claims, wherein the computer calculates the cardiac output contin- 
uously together with the heart rate HR from the stroke volume SV according to CO = HR • SV . 

24. The device as claimed in claim 23, wherein the computer determines the mean pressure MAP continuously from 
40 the blood pressure curve and thereby calculates the systemic resistance continuously according to 



45 



50 



0 MAP - CVP 
CO 

for an arbitrary central venous pressure CVP which has been measured or estimated. 

25. The device as claimed in claim 23, wherein the computer determines the mean pressure MAP continuously from 
the Wood pressure curve and thereby calculates the compliance function continuously according to 

C(p) MAP 2 CO 1 



( PW> 3 -MAP *p-3 *MAP 2 -p 



2 



55 
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(57) A device for individual in-vivo determination of 
the compliance function C(p) = dV/dp of the vascular 
system downstream of a ventricle and/or systemic blood 
flow of a living being from the blood pressure p(t) and a 
reference cardiac output COref. The device comprises 

a) a pressure sensor which continuously deter- 
mines the pressure p(t) in the aorta or in the vicinity 
of the aorta; 

b) a computer which 

b1) calculates the mean blood pressure MAP 
from the blood pressure p(t), 
b2) calculates the systemic resistance R of the 
body as 

MAP - CVP 
" COref ' 



CVP being an arbitrary central venous pres- 
sure which is ascertained or estimated, and 
COref being a reference value for the cardiac 
output, 

b3) takes at least the first differential of the 



blood pressure with respect to time 
p(t) = dp/dt , and 

b-4) calculates the compliance function C(p) at 
least from p(t), p(t) and R using a nonlinear 
model. 
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